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Abstract 

We have studied the chiral phase transition in three-dimensional QED in the 
presence of a Chern-Simons term for the gauge field. There exists a phase 
where the chiral symmetry is broken dynamically and we have determined the 
critical line for this symmetry breaking as a function of the effective coupling 
and the strength of the additional Chern-Simons term. In the presence of a 
Chern-Simons term, the chiral phase transition turns out to be first order, in 
sharp contrast to the phase transition in pure three-dimensional QED. 
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Dynamical symmetry breaking in three-dimensional quantum electrodynamics (QED3) 
has attracted much attention over the last ten years, both from a purely field theoretical 
point of view and because of its applications to condensed matter physics in connection with 
phenomena occurring in planar surfaces |]-[|]. A natural extension of pure QED3 is to add 
a Chern-Simons (CS) term for the gauge field ||, which breaks parity explicitly Indeed the 
statistics-changing CS term, together with the question whether or not there is a dynamically 
generated fermion mass, plays a key role for variants of QED3 to be effective theories for 
high-T c superconductivity and the fractional quantum Hall effect ||. Furthermore, QED3 
also has implications for high energy physics and physics of the early universe, since three- 
dimensional models are the high temperature limit of the corresponding four- dimensional 
theory. Recently, it has been suggested that the effective potential for high temperature 
QCD is also related to CS gauge theories f7|. 

The existence of the CS term leads to a novel feature in QED3, namely a first- order chiral 
phase-transition, as we show in this letter. The CS term generates a parity odd mass term for 
the fermions, but in addition there might be a parity even mass, which breaks chirality |||9] . 
Dynamical chiral symmetry breaking, a nonperturbative phenomenon, can be studied using 
the Dyson-Schwinger (DS) equation for the full fermion propagator. Both a numerical study 
of the full (truncated) DS equation and an analytical study of the approximated equations 
show that there is a first-order phase transition. This is quite remarkable and in sharp 
contrast to the infinite order phase transition one finds using the same truncation scheme 
in pure QED |,|. 

The Lagrangian in Euclidean space is 

£ = ip(i $ + e j/L- m e - Tm )ip + \F 2 V 

+ 2^ 6 m^p^^^p + ^gauge fixing > W 

with the dimensionful parameter 9 determining the relative strength of the CS term. We 
use four-component spinors for the fermions, and a four- dimensional representation for the 
7-matrices. The matrix r is defined in such a way that the term m %pTip is odd under a 
parity transformation [p],[T0|| . Also the CS term is odd under a parity transformation, the 



other terms in the Lagrangian are invariant under a parity transformation. 

With such a representation we can define chirality similar as in four-dimensional QED. 
Without an explicit mass m e for the fermions, the Lagrangian is chirally symmetric, but 
the parity even mass m e breaks this symmetry. Note that the other mass, m , is chirally 
invariant. Just as in pure QED, chiral symmetry can be broken dynamically due to nonper- 
turbative effects, which can be studied by solving the DS equation for the fermion propagator 
with both explicit masses m e and m equal to zero. 



The standard way to truncate the DS equation in QED3 is the 1/N expansion [11], where 
iV is the number of fermion flavors. The coupling constant e 2 has the dimension of mass, 
and we use the large N limit in such a way that e 2 j and the product N e 2 remains fixed: 
Ne 2 = 8a with a fixed. In this 1/N expansion the one- loop vacuum polarization has to be 
taken into account, because this vacuum polarization is of order one: using bare massless 
fermions, the transverse part of the vacuum polarization is just n T (g) = — a\q\ [|l|. It is easy 
to show that at order one there is no parity odd part of the vacuum polarization. 

The full vertex is replaced by the bare one, because that is the leading order contribu- 
tion in 1/N. In order to be consistent with the requirement that the vertex renormalization 
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and the fermion wave function renormalization are equal, we use a suitable nonlocal gauge 
function In pure QED, one can construct a gauge in which the wave function renormal- 



ization is exactly one. In the presence of a CS term, this condition can only be satisfied up 
to order 9/N, but we are considering small 9 only. The proper choice for the gauge function 
is i 

( 9 n 2a 4a9 9\q\ 
a(q) = 2 q D (<?) + -j— r H — arctan- 



\q\ q 2 a 2 + a\q\ + 9 2 

a 2 -9 2 , a 2 + 9 2 \ 
+ ^ 2 - ln (a + \q\) 2 + 9 2 - 1 )- (2) 

With this gauge, we also satisfy the Ward-Takahashi identity up to corrections proportional 
to 9/N and to the dynamically generated mass function, which are both negligible. Gauge 
covariance can (in principle) be recovered by applying the Landau-Khalatnikov transforma- 
tion rules to the various Green's functions 

The inverse full fermion propagator can be written as 

S^ip) = A e (p) p , + A (p)r 0-B e (p) - B (p)r. (3) 

The functions A(p) and B(p) are scalar functions of the absolute values of the momenta, and 
their bare values are A e — 1, A — 0, B e = m e , and B Q = m . We use the decomposition 
A± = A e ± A a and B± = B e ± B Q , which leads, together with the above truncation scheme, 
to the following two sets of coupled integral equations 

A ■ to f d3fc 2B ± (k)D°(q) p-q 

±{P) Np 2 J (2tt) 3 k 2 A 2 ± (k) + B 2 ± {k) \q\ ' K ' 

d 3 k 1 

■x 




(2^fk 2 A 2 ± {k) + B 2 ± {k) 

^, a ^)^2A ±{ k)DO {q) k -f 
\ Q J 111 



h(k) [ 2D 1 (q) + ^) T 2A ± (k)D°(q)^ ) , (5) 



where D T and D° are the transverse and the parity odd part of the gauge boson propagator 



(q 2 + a\q\) 2 + 9 2 q 
= (q 2 + a\q\) 2 + 9 2 q 2 ' 

and q = k — p. Note that the equations for A + and B + decouple from the ones for A_ and 
B_. It is also important to observe that once we have found a solution for A + and B + , we 
automatically have also a solution for A_ and B_: namely the set A_ = A + and B_ = —B + . 
That means that we can always construct a chirally symmetric (but parity odd) solution, 
with B e = 0. The question of dynamical chiral symmetry breaking turns into the question 
whether or not there exist two (or more) solutions of the set of integral equations. 

Without the CS term there is dynamical chiral symmetry breaking only for A > A c = 3/16 
||, where we have defined the effective coupling A = 8/(A% 2 ). We expect a similar situation 
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in the presence of the CS term, at least if the parameter 6 is small. That means that for 
A < A c we only have the chirally symmetric solution of the above equations, but for A > A c 
we expect that there are (at least) two solutions for both B + and B_ possible, in such a way 
that there is a nonzero solution for B e . An essential difference from pure QED is that in the 
presence of the CS term there is no trivial solution B±(p) = 0. Due to the explicit breaking 
of parity, the fermions always acquire a parity-odd mass term B Q , even if the explicit odd 
mass term m is zero. 

Firstly, we solve the DS equation analytically after some further approximations. Using 
A(jp) = 1 + O{0) due to the nonlocal gauge, we replace A{p) by one, so we get an integral 
equation for B± only, consisting of two terms, Eq. (|5]) with A±(jp) = 1. The first term is the 
same as in pure QED3, and the essential region for this term is the infrared fl]], p, k << a. 
So we consider the integral for small momenta only, expand the integration kernel in powers 
of p and k, and introduce a cutoff at k = a. We also linearize the equation, by replacing the 
denominator k 2 + B±(k) by k 2 + B 2 _(0), which is reliable as long as B±(p) is almost constant 
for small momenta. In pure QED these approximations lead to almost the same result as 
the full nonlinear integral equation [JTjynj. 

The second term, proportional to 6, can be calculated by neglecting B± with respect to 
k 2 in the denominator, and expanding the integrand in powers of min(p, q) / max(p, q) and 
9. Taking into account only the leading order terms gives in the infrared region 

8a f d 3 k 2D°(q) k ■ q 



(2tt) 3 k 2 + Bl(k) \q\ 
±\9 + 0{p) + <D(6 3 ) , (8) 



and in the ultraviolet region p > a 

I 6 (p)^± 1 -^ + 0(l/p 2 ) + 0(6 3 ). (9) 

This means that in the ultraviolet region the CS term will dominate, since without the CS 
term, B(p) falls off much more rapidly in the far ultraviolet. Higher order contributions in 
min(p, q)/ max(p, q) will slightly change this result, but not affect the general behavior [ITS . 
Thus we have for p < a and to order 8 

B± {p) = f A f dk - fc2 _ B± ^\, ± Afl , (10) 
' k 2 + M| max(p, k) 



where we have defined M± = B±(0). This integral equation can easily be solved by con- 
verting it to a second-order differential equation with boundary conditions. The solution 
is 

B ± (p) = M± 2 F 1 {a+, a-, §; - P 2 /M 2 ± ) , (11) 



where a± = ^(1 ± zy'16A/3 — 1). The ultraviolet boundary condition leads to the condition 

M± 2 F 1 (a+,a-,l;-a?/Ml)=±\6. (12) 

In order to determine the behavior of M as a function of 6, we can plot the LHS of 
Eq. (|TJ) divided by A for a given value of A as function of the mass parameter M, see Fig.[l] 
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(for convenience we have set a = 1 in our figures, which just defines the energy scale). 
From this figure we can see that there are three solutions possible for B + and B_ at small 
values of 9. A closer look at the region around the origin would reveal that there exist more 
solutions for extremely small values of 9. In the absence of a CS term there are infinitely 
many oscillating solutions [I], but it has been shown that the vacuum corresponds to the 
nodeless one, with the highest value for \M\. With the CS term, there is only a finite number 



of oscillating solutions P,|T5 

The chirally symmetric solution consist of the combination of B + (p) and B_(p) with 
B_(p) = —B + (p), which can always be constructed. A solution which breaks chiral sym- 
metry can only be constructed if there is a different solution M_ of Eq. ([12]), corresponding 
to -B-(p) which is not equal to — B + (p). As we can see from Fig.|l|, this is only possible for 
9 < 9 C (X), beyond this critical value there is only one solution possible for B + and B_ which 
automatically gives B e = and B Q = B + = —B_. At the critical value M± does not become 
zero, nor does B e (0) = (5 + (0) + B_(0))/2, which can be regarded as the order parameter 
of the chiral phase transition. This clearly signals a first-order phase transition, in sharp 
contrast to the pure QED case. 

We can also plot M± versus A for a fixed value of 9, see Fig. |[ Here we see that if we 
increase A for fixed 9, the chiral symmetry breaking solutions appears if A exceeds some 
critical value X c (9), which increases rapidly as a function of 9. This figure shows that the 
chiral phase transition is first order in this direction as well: increasing A beyond A c ((9) gives 
rise to the second (and third) solution, but at the phase transition neither M± nor B e (0) 
become zero. In this figure we can also see that in the limit 9 — > both M± and M± 
go towards the nontrivial solution m of pure QED, and the critical value A c goes towards 
X c (9 = 0) = 3/16. 

The critical parameters A c and 9 C can be calculated directly from Eq. |T2[ In Fig. ^] we 
have shown the critical line in the (A, 6 l )-plane. For small values of M we can make an 



expansion in order to get an explicit expression for 9 C (X). To leading order in ^/l6A/3 
this gives 



9 ~exp(-37r/^16A/3- 1) . (13) 

Secondly, we have solved numerically the set of coupled integral equations for A and B, 
Eqs. |] and |5|, without any further approximations, and these numerical results are qualita- 
tively in good agreement with our analytical results. First we checked our assumption that 
A±(p) is close to one, and it turns out that the deviation is indeed negligible for small values 
of 9 ]T5|] . Furthermore, we have found numerically the following solutions for B + , using the 
notation m(p) for the solution with 9 = 0: 

1. B+(p) = 0(m(p)) > for A > X c (9 = 0), 
B + {p) = 0{9) for A < X c {9 = 0); 

2. B+(0) = C(-m(0)) < 0. 

The first solution exists for all values of both A and 9, whereas the second one exists only 
for values of A > A c (#) and 9 < 9 C (X). Note that this second solution behaves like 9/p in 
the ultraviolet, so it has some node at a particular value of p. The existence of this second 
solution allows for a nonzero chiral symmetry breaking solution B e (p) = (B + (p) + B_(p))/2. 
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Numerically, it is extremely difficult to establish a first-order phase transition and to 
determine the critical values A c and 9 C . However, our numerical results all support our 
analytical results, and indicate strongly that the chiral phase transition is indeed first order. 
In Fig. |] we have shown the behavior of B + (0) and B_(0) at fixed 9 as a function of A. We 
can see that the behavior is the same as in Fig. |2|: increasing A at fixed 9 leads to a second 
solution B± beyond some critical value A c > X c (9 = 0). Close to the critical value B±(0) 
does not go to zero, nor does B e (0), signaling a first-order chiral phase transition. Also the 
behavior for increasing 9 at fixed A is qualitatively the same as our analytical result. 

In conclusion, both the numerical and the analytical results show that there is a first-order 
chiral phase transition in QED3 with explicit CS term. This result is very remarkable, given 
the well-known infinite order phase transition (Miransky-scaling) |l| in the absence of the 
CS term. Also the other known chiral phase transitions in four-dimensional gauge theories 
are of second (or higher) order. This first-order phase transition is a new and interesting 
phenomenon, and it might lead to new insights into chiral phase transitions in general. In 
particular, the connection between the CS term and the first-order phase transition should 
be studied in more detail. 

This result should also be contrasted with some previous results in analyzing this model 
||, indicating just a minor quantitative effect on the critical coupling and scaling behavior 
due to the CS term. Both our numerical and analytical results reveal that the presence of 
an explicit CS term changes the nature of the chiral phase transition drastically. 
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and he would like to thank the members of the Graduate School of Science and Technology 
for their hospitality during that stay. We would like to thank K. Yamawaki, Yoonbai Kim, 
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JSPS (P. M. being a JSPS fellow under number 94146). 
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FIGURES 



FIG. 1. 6 as function of M for some different values of A. The upper half of the plane 
corresponds to solutions for B + , the lower half for £>_. Note that both the upper and the lower 
half correspond to positive values of 6. 

FIG. 2. \M±\ as function of A for some different values of 6. 

FIG. 3. The critical line for the chiral phase transition in the (A, #)-plane. 

FIG. 4. The infrared values -B+(0) and -B_(0) as function of A for some different values of 6. 
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